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Abstract — We give an overview of tlie work done witli the 
LNE-SYRTE fountain ensemble during tlie last five years. After 
a description of the clock ensemble, comprising three fountains 
FOl, F02 and FOM, and its newest developments, we review 
recent studies of several systematic frequency shifts. This includes 
the distributed cavity phase shift which we evaluate for the 
FOl and FOM fountains, applying the techniques of our recent 
work on F02. We also report calculations of the microwave 
lensing frequency shift for the three fountains, review the 
status of the blackbody radiation shift, and summarize recent 
experimental work to control microwave leakage and spurious 
phase perturbations. We give current accuracy budgets. We also 
describe several applications in time and frequency metrology: 
fountain comparisons, calibrations of the international atomic 
time, secondary representation of the SI second based on the 
*^Rb hyperflne frequency, absolute measurements of optical 
frequencies, tests of the T2L2 satellite laser link, and review 
fundamental physics applications of the LNE-SYRTE fountain 
ensemble. Finally, we give a summary of the tests of the PHARAO 
cold atom space clock performed using the FOM transportable 
fountain. 



I. Introduction 

Atomic fountain clocks provide the most accurate realiza- 
tion of the SI second and define the accuracy of the widely 
used international atomic time (TAI). Ten years ago, a few 
formal TAI calibration reports were available each year. Now, 
there are normally several formal reports each month as a 
result of a great deal of effort in several metrology institutes 
to improve the reliability and long term operation. At LNE- 
SYRTE, the fountain ensemble typically provides more than 
20 TAI calibrations each year, about half of all reports. These 
large improvements in reliability and corresponding capability 
to perform lengthy measurements yield an enhanced ability to 
investigate smaller and smaller frequency shifts, to compare 
fountain clocks with more stringent uncertainties and to test 
their long term behavior In turn, applications using atomic 

Manuscript received September 15, 2011; accepted January 23, 2012. 

This work was supported by SYstemes de Reference Temps-Espace 
(SYRTE), the Laboratoire National de Metrologie et dEssais (LNE), the 
Centre National de la Recherche Scientifique (CNRS), Universite Pierre et 
Marie curie (UPMC), the Observatoire de Paris, the Institut Francilien de 
Recherche sur les Atomes Froids (IFRAF), the National Science Foundation 
(NSF), Penn state, the Ville de Paris, the European Space Agency (ESA), and 
the Centre National dEtudes Spatiales (CNES). 

J. Guena, M. Abgrall, D. Rovera, P. Laurent, B. Chupin, M. Lours, G. 
Santarelli, P. Rosenbusch, A. Clairon and S. Bize are affiliated to: LNE- 
SYRTE, Observatoire de Paris, CNRS, UPMC, 61 avenue de I'Observatoire, 
75014 Paris, France 

M. E. Tobar is affiliated to: The School of Physics, University of Western 
Australia, Crawley, Western Australia 

R. Li and K. Gibble are affiliated to: The Department of Physics, The 
Pennsylvania State University, University Park, Pennsylvania 16802, USA 

This work was supported by SYRTE, LNE, CNRS, UPMC, Observatoire 
de Paris, IFRAF, NSF, Penn State, Ville de Paris, ESA and CNES 

DOT: [http://dx.doi.org/10. 1 109/TUFFC.2012.2208| 



fountain clocks have also developed and benefited from the 
improvements of the last few years. 

In this article, we give an overview of developments and 
of applications of the LNE-SYRTE atomic fountain ensemble 
during the last five years. 

II. LNE-SYRTE FOUNTAIN ENSEMBLE 
A. Overview 

LNE-SYRTE has operated three atomic fountains for more 
than a decade. The first, FOl, is a ^^'^Cs fountain which has 
been in operation since 1994 fT\. The second, FOM, was 
originally a prototype for the PHARACQ cold atom space 
clock 121 and was later modified to be a transportable fountain 
clock. The third, F02, is a dual fountain which operates with 
^^Rb and ^^'^Cs simultaneously [3J. Development for more 
than a decade has continuously improved each of these. Our 
description here reflects the present configuration. 

The three fountains share several features. All of them are 
gathering atoms in a Lin _L Lin optical molasses with laser 
beams oriented to launch atoms vertically in the so-called 
(1,1,1) configuratiorj^ |4|. The 6 laser beams for the optical 
molasses are obtained from 6 fiber-coupled collimators, which 
are designed, machined and tuned to ensure proper alignment 
at the few 100 /irad level when placed against the correspond- 
ing reference surfaces of the vacuum chamber Similar fiber- 
coupled collimators provide the other laser beams for state 
selection and detection. In all three fountains, including the 
^^Rb part of F02, the laser system is on a separate optical 
bench and the light is coupled to the vacuum system with 
optical fibers. The lasers are home-built extended cavity diode 
lasers, most of which use a narrow band interference filter for 
wavelength selection f5l. This design is used in the PHARAO 



cold atom space clock (see VI i. Also, it has recently become 
commercially available. The required power is obtained either 
by injection locking another laser diode or using a tapered 
semi-conductor laser amplifier In all three fountains, atoms are 
launched with the moving molasses technique [61 at a velocity 
of ^ 4 m.s~^ and a temperature of ^ 1 /iK. Atoms exit the 
launch sequence in the upper hyperfine state. For normal clock 
operation, microwave and laser pulses further select the atoms 
in the mp = Zeeman sub-state of the lower hyperfine state 
(\F = 3, = 0) for "^Cs, |F = 1, wj. = 0) for s'^Rb). 

All three fountains also have at least one layer of magnetic 
shield that surrounds the entire vacuum system. Notably, the 

'PHARAO stands for "Projet d'Horloge Atomique par Refroidissement 
d' Atomes en Orbite" 

^In this configuration, the 3 pairs of counter-propagating laser beams for 
the molasses are aligned along the axes of a 3 dimensional orthonormal basis, 
where the (111) direction is along the vertical direction, as sketched in Fig.|2] 
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molasses and the detection region are shielded, which is 
important to operate in the presence of the strong magnetic 
field fluctuations as of many urban environments, including 
the Observatoire de Paris. The largest, vertical magnetic field 
component in the lower part of the fountain (molasses and 
detection) is measured by a flux-gate magnetometer and ac- 
tively stabilized with a set of horizontal coils distributed over 
the height of the fountain. Two or three additional cylindrical 
shields with end caps surround only the interrogation region 
to further attenuate magnetic field fluctuations. Inside the 
innermost shield, a solenoid and a set of compensation coils 
provide a static magnetic field in the range of 80 to 200 nT. 
The field is homogeneous to 10^"^ and stable to 2 x 10^^ 
from 1 s to several days, as established by spectroscopy of the 
\F,mF = 1) — > \F + l,mF = 1) field sensitive transition. 
This corresponds to a stability of ~ 4 pT. 

All three fountains also use the same strategy to control 
temperature and thermal gradients in the interrogation region. 
The innermost layer surrounding the interrogation is made of a 
material with high thermal conductivity (copper or aluminium 
alloy) and is well isolated from the envkonment to ensure good 
temperature uniformity of this layer. The temperature of this 
innermost layer is monitored using platinum resistors and is 
allowed to drift slowly, following temperature fluctuations in 
the room with a typical time constant of 1 day. The measured 
temperature is used to infer the blackbody radiation tempera- 
ture seen by the atoms. With this scheme, the temperature of 
the innermost layer is uniform to much better than 0.1 K and 
no large gradients exist between the inner region and the envi- 
ronment. This essentially removes the need to care about the 
impact of the necessary openings (for the atomic trajectories, 
pumping,...) on the effective blackbody temperature for the 
atoms. Using heating devices to change or actively stabilize 
the temperature of the innermost layer induced a significant 
degradation of the temperature uniformity, which was difficult 
to compensate, even with segmented heaters. Consequently, 
this approach was abandoned. 

The microwave interrogation cavities in all three fountains 
are TEqh copper resonators with two independent microwave 
feedthroughs. The microwave synthesizers used to feed the 
cavities are home built. Despite significant differences, these 
synthesizers share the capability to use the ultra-low phase 
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noise reference signal (see below II-Ei from a cryogenic sap- 
phire oscillator (CSO) without significant degradation at the 
level of one to a few parts in 10^^ at 1 s [iTl. Consequently, the 
short-term instability of all three fountains is typically limited 
by the atomic quantum projection noise [8| and technical noise 
in the atom detection. The best short-term stability to date for 
the operation of an atomic fountain with a high atom number 
is 1.6 X 10~^'^T~^/'^ as reported in |9| and exemplified in Fig. 
[T] The present version of the synthesizers also incorporate a 
phase-stable microwave switch |10| which helps to assess and 
remove leakage fields that may affect the interrogation process, 
independent of other systematic effects. Some details are given 
in HITDl 

For all three fountains, the computer system allows auto- 
mated control of most of the important parameters, enabling 
routine implementation of multiple interleaved measurements 
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Fig. 1. Fractional frequency instability of F02-Cs fountain when using 
a cryogenic sapphire oscillator as a reference signal. The black points are 
obtained for a high atom number: Nat = 5 X 10® detected atoms. The red 
points are obtained for Nat/'i detected atoms. The blue line corresponds to 
1.6 X 10~^*T~^/^. Above 100 s, frequency fluctuations of the ultra-stable 
interrogation signal (see |II-E| , common to the high and low atom number 
measurements, are apparent. 



with different configurations to study systematic shifts. A 
large number of parameters (laser powers, temperatures in the 
interrogation region, temperatures in the environment, at the 
Cs or Rb sources, humidity,...) are frequently monitored. The 
microwave power is regularly optimized and the magnetic field 
is verified via spectroscopy of a field sensitive transition, with 
automated computer controlled sequences. 

B. FOl 

FOl uses a 2 dimensional magneto-optical trap (2DMOT) 
ifTTI as a source of slow atoms to load the optical molasses. 
An advantage of this type of source is that it confines the 
high alkali vapor pressure to the 2DMOT chamber, yielding 
a large number of cold atoms in the molasses while keeping 
source consumption low. It also facilitates a better vacuum in 
the molasses region, thereby minimizing perturbations in the 
detection region and the contamination of the interrogation, 
including the microwave cavity, by background alkali vapor 
With higher laser powers in the 2DMOT, it should be possible 
to reproduce the atom numbers previously obtained with 
a chirped-slowed thermal atomic beam, that gave the best 
reported fountain stability |9|, while consuming an order of 
magnitude less of alkali metal. So far, the laser power devoted 
to the 2DMOT yields ^ 4 times fewer atoms in the molasses 
than with our best chirped-cooled atomic beam. Note that 
2DMOT loading of the molasses has some disadvantages. 
Besides the quite intricated mechanical and optical design and 
the need for extra laser sources, we observe a degradation 
of the 2DMOT laser windows that are in contact with the 
relatively high alkali vapor pressure (up to 10~^ Pa) after a 
few years of almost continuous operation. Also, loading the 
molasses asymmetrically from the slow atom beam from the 
2DMOT, leads to a non-spherical, complicated atomic cloud 
distribution which makes the investigation of some systematic 
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frequency shifts more difficult. Note however, that by using the 
adiabatic passage method to change the atom number lfT2l . 
the extrapolation of the cold collision shift to zero density 
is essentially unaffected by distortions of the atomic cloud's 
shape. 

The microwave cavity in FOl has a diameter of 42.99 mm 
and a length of 42.88 mm. The below-cutoff guides that let 
the atoms pass through the cavity have a diameter of 12 mm 
and a length of 60 mm. At the two ends of the cutoff guides, 
an additional diaphragm with a diameter of 11 mm blocks 
atoms that could travel close to the copper walls. Each of 
the 2 opposing microwave feedthroughs includes a resonant, 
polarization-filtering rectangular waveguide section. In the 
latest implementation, the cavity has a loaded quality factor 
of 14000. 

The interrogation region in FOl is currently equipped with 
Stark plates. These can be used to measure the sensitivity of 
the clock transition to static electric fields to determine the 
Stark coefficient for the blackbody radiation shift correction. 
Alternatively, FOl can be equipped with a blackbody radiator 



to directly measure the blackbody radiation shift (see III-C 
below). 

C. FOl dual Rb/Cs fountain 

A schematic view of the F02 dual Rb/Cs fountain is shown 
in Fig. |2] It uses 2 independent 2DM0T sources for Cs 
and for Rb. For Rb, only the 27.8% abundant isotope is 
used in the fountain and the use of a 2DMOT eliminates 
the ^^Rb background atoms, as for Cs, and also the larger 
background from the unwanted *^Rb isotope. The Rb and Cs 
molasses are superimposed with dichroic collimators for each 
molasses beam, where the 852 nm light (Cs) and the 780 nm 
light (Rb) come from 2 independent optical benches through 
polarizing fibers and are then combined with a dielectric 
dichroic beamsplitter and collimated with an achromatic lens 
lfT3l . Similar dichroic collimators are used for all other laser 
beams, the pushing beam for the state selection, detection 
beams and repumper for detection. 

As shown in Fig. [2j there are a total of 4 microwave cavities 
for state selection and interrogation of Rb and Cs atoms. 
Because both interrogation cavities are in the interrogation 
region, their temperature must be uniform and therefore the 
same as the environment temperature, which removes any 
possibility to independently tune each cavity. Practically, the 
Rb and Cs cavities are a single copper assembly where the Rb 
and Cs cavity resonances were mechanically tuned relative to 
one another, after a lengthy and tedious process. The resonance 
frequencies had to be verified under vacuum and tightly 
controlled thermal conditions, after each tuning step that was 
a controlled polishing of the two cavity spacers |14||l31. The 
Rb and Cs microwave cavities have the same aspect ratio, 
with the radii scaled by the ratio of the hyperfine frequencies. 
The cavity for Cs is located at the top. It has a diameter of 
50.00 mm and a length of 26.87 mm. The Rb cavity has a 
diameter of 67.25 mm and a length of 36.01 mm. The below- 
cutoff end cap tubes have a diameter of 12 mm and a length of 
~ 50 mm. The 2 opposing microwave feedthroughs for each 



cavity include a resonant, polarization filtering rectangular 
waveguide section. The loaded quality factor is 7000 for Cs 
and 6000 for Rb. 

The Rb and Cs clocks in F02 routinely operate simulta- 
neously since 2009 |3|. Two independent but connected and 
synchronized computer systems allow the coordination of mul- 
tiple interleaved measurements with Cs and Rb, as well as all 
other tasks for each clock, such as monitoring the temperature, 
the magnetic field, etc. Thus far, the nominal configuration 
is to launch the Rb and Cs clouds almost simultaneously, 
but at different velocities, so that the interrogation times are 
simultaneous, i.e. the Rb cloud is at the center of the Rb cavity 
when the Cs cloud is at the center of the Cs cavity, on both 
the way up and the way down. In this way, the two atomic 
clouds do not interact with each other during the interrogation 
process, avoiding interspecies collisions. The two clouds also 
remain well separated as they fall to the detection region, 
allowing the time-resolved detection of Rb and Cs |3|. As 
a result, the Rb and Cs clocks run simultaneously without 
impacting the performance of the other. The flexible computer 
systems allow for a large number of other possibilities which 
will be explored in the future. For instance, one cloud can be 
purposely launched through the other with a well controlled 
and yet variable speed, in order to study interspecies collisions 
at varying energies llT6l . based either on measurements of the 
cold collision frequency shift or on direct detection of the 
scattered wave I.17j. 
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Fig. 2. Schematic view of the F02 dual Rb/Cs fountain. The 852 nm light 
and the 780 nm light for ^'^'^Cs and **^Rb atoms are superimposed in each 
of the 9 laser beams with dichroic beam splitters. ^'^^Cs and *^Rb atoms are 
captured, launched, state selected, probed and detected simultaneously at each 
fountain cycle of 1.6 s. There are 2 independent 2D magneto-optic trap for 
the two atomic species. 
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D. FOM transportable fountain 

Photographs of the transportable fountain FOM are shown 
in Fig. [3] FOM has two major subsystems. The first is the 
vacuum system where atoms are manipulated. The second 
is the optical laser bench, its electronics, the microwave 
synthesizer and the computer system that operates the fountain, 
all of which are mounted on a single frame. Both of these 
subsystems include vibration damping devices necessary from 
transportation. When disconnected, the two subsystems fit 
into a small size truck. A battery powers the ion pump 
to maintain the vacuum of the system. After transportation, 
the two subsystems are reconnected (polarizing optical fibers 
for the laser light, microwave cables, wires for monitoring) 
and the system is allowed to thermally stabilize for a few 
hours. After a limited number of manual re-optimizations 
(laser power through the optical fibers), a computer controlled 
sequence is started which automatically finds the laser set 
points and activates the laser stabilization loops fTS^. At 
remote locations, close-to-nominal operation of the fountain 
is typically recovered in less than 2 days. 

One notable feature of the transportable fountain FOM 
is that the microwave synthesizer has two possible modes 
of operation. At LNE-SYRTE and at laboratories where a 
high quality metrological reference signal is available, the 
microwave synthesizer is synchronized to this external refer- 
ence and the FOM fountain works as a frequency standard 
calibrating the frequency of the external reference. In the 
second configuration, FOM autonomously delivers a primary 
reference signal. Here, a BVA quartz oscillator synchronizes 
the microwave synthesizer and a digital loop locks the fre- 
quency of this BVA quartz oscillator to the atoms. 

The microwave interrogation cavity in FOM has a diameter 
of 41.06 mm and a length of 65.66 mm. The two opposing 
microwave feedthroughs have direct evanescent coupling: the 
coaxial input cable is terminated in a small non-resonant 
cavity, with an antenna located a few millimeters away from 
the coupling iris of the cavity. The cavity has a loaded quality 
factor of 21000. 

E. Ultra-stable reference and frequency distribution scheme 

Figure |4] is a schematic of the LNE-SYRTE fountain en- 
semble. A Cryogenic Sapphire Oscillator (CSO) |fT9lll20l is 
the local oscillator for all of the atomic fountains. Its ultra- 
low phase noise allows the fountains to operate at the quantum 
projection noise Umit lH), producing short-term instabilities as 
low as 1.6 X 10~^^ at 1 s (f9l and Fig.[l]i. As shown in Fig.|4] 
a finely tunabl^ microwave 11.98 GHz signal is first derived 
from the free-running CSO. This frequency is divided down 
to 100 MHz and compared to the output of one of the LNE- 
SYRTE hydrogen masers, HI in Fig. |4] A computer driven 
digital phase-lock loop (time constant ^ 1000 s) corrects 
the frequency of the 11.98 GHz microwave signal to phase 
lock the 100 MHz signals. This 11.98 GHz reference has the 
exquisite phase noise properties of the CSO for time scales 
less than ^ 1000 s and the mid and long-term frequency 

'with a computer controlled Direct Digital Synthesizer 



Fig. 3. Pictures of the transportable fountain FOM. On the left is the vacuum 
chamber for cooling and launching cesium atoms. Four layers of magnetic 
shielding surround the chamber. On the right is the second clock component 
with on top the optical bench inside an acoustic box, and below the electronics, 
including the microwave synthesis chain. A computer controls the clock. 



Stability of the hydrogen maser, while providing the connection 
to LNE-SYRTE timescales, to the International Atomic Time 
(TAI) and to other laboratories via GPS and TWSTFT satellite 
time transfer This 11.98 GHz signal is the high-performance 
reference from which all other signals are derived. A 1 GHz 
source is the starting point of the microwave synthesis for 
FOl, F02-Rb and FOM. Since F02 is located next door 
to the CSO, a more direct synthesis from 11.98 GHz to 
9.192 GHz could be used for F02-Cs ||2T1 . To compare to 
optical clocks through optical frequency combs, a 8.985 GHz 
is also generated. Finally, an additional synthesizer generates 
signals in the 6 to 9 GHz range for testing (see |III-D| i. The 
design and characterization of the low phase noise electronics 
that generates these reference signals are reported in |21||7j. 
Phase noise power spectral densities can be found in II20I 
for the CSO alone, and in 121J, ITJ and ||22J for the derived 
reference signals. 

Figure |4] shows 100 to 300 m long actively compensated 
optical fibers links that are used to distribute the reference 
signals. These links are simplified versions of similar long- 
distance links 1231 . due to the short distance and limited 
insertion losses. The reference signal modulates the current 
of a 1.55 jjm diode laser that is launched into a standard 
telecom fiber. At the remote end, a fast photodiode detects the 
intensity modulation and delivers the signal to the application. 
A fraction of this recovered signal is used to modulate a second 
diode laser that is launched back into the fiber At the emission 
site, the returned signal is detected by another fast photodiode 
and then mixed with the reference signal to detect the phase 
variations imposed by the two passes propagation through the 
fiber A temperature controlled fiber spool of 50 to 100 m 
is inserted to control the length of the link to cancel these 
phase variations with a computer operated servo-loop (time 
constant ~ 10 s). For the high performance links at 6 — 9 GHz 
and 8.985 GHz, an additional fast actuated piezo-driven fiber 
stretcher is used (bandwidth ~ 400 Hz). Characterization of 
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these links has shown stabiHties ^ 10^^'* at 1 s for the 1 GHz 
caiTier and ~ 2 x 10~^^ at 1 s for the 6 — 9 GHz and 
8.985 GHz carriers, and a residual frequency offset well below 
10^^^. A computer acquires and analyzes the doubled-passed 
phase signals to assess the performance in real-time and detect 
anomalies. 
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Fig. 4. LNE-SYRTE clock ensemble. A Cryogenic Sapphire Oscillator 
generates an ultra-low noise reference signal that is distributed to all the three 
fountains. This 11.98 GHz reference signal is slowly phase-locked (PLL) to 
the output of a hydrogen maser (HI), which provides long term stability 
and the connection to local timescales, satellite time transfer systems (GPS, 
TWSTFT), and the international atomic time (TAI). The reference signal can 
also be distributed through actively compensated fiber links. 



F. New solution for the ultra-low noise reference oscillator 

The CSO allows a large improvement in the fountain short- 
term stability but requires cumbersome and costly liquid 
helium refills every 26 days. Each refill lasts several hours and 
degrades the performance of the CSO for nearly half a day. An 
alternate solution for fully continuous operation of an ultra- 
low noise microwave reference would be to use a pulsed-tube 
cryocooler for the CSO |24||25||26||27|. However, advances 
in laser frequency stabilization and optical frequency combs 
offer the possibility to replace the CSO with a fully continuous, 
cryogen-free and maintenance-free ultra-low noise microwave 
generation system. 

Figure|5]shows a schematic of such a laser-based system. An 
ultra-stable laser is realized by stabilizing a reliable 1.55 jim 
laser, such as an Erbium-doped fiber laser or an extended 
cavity laser, to an ultra-stable reference Fabry-Perot cavity. 
Several Fabry-Perot stabilized lasers have demonstrated short- 
term stabilities well below 10"^^ ll28lll29'H30|, similar to 
stabilities reported for cryogenic sapphire oscillators ll20l . 
Note however that simplified versions of these cavities can 
be sufficient for application to atomic fountains. Microwave 
signals with short-term frequency stabilities in the mid 10~^^ 
range at 1 s are indeed sufficient to achieve the quantum 



projection noise limit at 10^^^ at 1 s. This instability is not 
easy to surpass for other reasons: it is difficulty to increase 
the number of useful atoms per cycle, to cope with the 
tremendous cold collision shift that would result from more 
atoms, and to suppress certain detection noise sources. A short 
(< 10 cm) reference Fabry-Perot cavity with mirrors made 
of Ultra Low Expansion glass with a comparatively simple 
temperature stabilization system can deliver instabilities less 
than 2 x 10^^^ between 1 and 1000 s, which is well suited to 
atomic fountains. 

As shown in Fig. |5] the ultra-stable laser is used to stabilize 
the repetition rate of the optical frequency comb generated by 
a fiber based femtosecond laser. A well-chosen harmonic, here 
12 GHz, is detected and divided down to 100 MHz to allow a 
phase comparison with a hydrogen maser, similar to what is 
done for the CSO (see |II-E| l. Here, a digital phase lock loop 
shifts the frequency of the acousto-optic modulator (AOM) and 
thereby the ultra-stable light that controls the comb repetition 
rate, ensuring the long-term stability of the 12 GHz reference 
that ultimately provides the 9.192 GHz interrogation signal for 
the fountains. 

Ultra-low noise microwave generation with fiber-based op- 
tical frequency combs was demonstrated in f3Tl and fy2\. 
Atomic fountains using microwave signals derived from op- 
tical references reported fountain short-term stabilities of 
3.5 X 10-i4r-i/2 |33] and 7.4 x IQ-^^t'^/^ For both, 
the fountain stability was quantum projection noise limited and 
the microwave signal alone was compatible with a fountain 
stability < 10~^'*t~^/^. Subsequent progress on fiber-based 
optical frequency combs will enable the implementation of a 
fully operational version of the system of Fig. |5] 

Note that an alternative approach for achieving quantum 
noise limited stability without an ultra-low noise microwave 
sources is the continuous atomic fountain If35lll36l . 

III. Recent studies of systematic shifts 

A. Distributed cavity phase shift 

The distributed cavity phase (DCP) shift is a residual first 
order Doppler shift that occurs when the moving atoms interact 
with the field inside the microwave cavity that has a spatial 
phase variation. Up until recently, measurements have been 
unable to reveal this effect and to reasonably agree with cal- 
culations. On the one hand, this was because the corresponding 
frequency shifts can be quite small, the measured frequency 
shifts give only indirect information about the phase variations, 
this shift and its changes with fountain parameters can be 
intricated with other systematic effects. On the other hand, 
accurately calculating the actual phase distribution and its 
corresponding frequency shift was a challenging numerical 
problem. 

In ||37]| and Il38l . a new theoretical approach was proposed to 
solve the problem of the cavity field computation. The field, 
represented explicitly as the sum of a large standing wave 
Ho(r) and of a small field, g(r) ^ g,„(p, z) cos(m(/)), 
that describes the effects of wall losses and cavity feeds, is 
expressed as an azimuthal Fourier series. Only the lowest 
azimuthal terms are relevant m = 0, m = 1 and m — 
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Fig. 5. An ultra-stable laser and an optical frequency comb produce an 
ultra-low noise microwave flywheel for the atomic fountains. The 12 GHz 
microwave signal is obtained by detecting a harmonic of the repetition rate of 
a Er-doped fiber femtosecond laser, which is stabilized to a CW laser itself 
stabilized to an ultra-stable Fabry-Perot cavity. 



2. We have recently published measurements of the DCP 
shift performed with ^■^^Cs in F02 1391 where, for the first 
time, quantitative agreement between theory and experiment 
is found. This study led to a better characterization and a 
quantitative evaluation of the contributions to the uncertainty 
of the DCP shift. It highlighted the importance of controlling 
the fountain geometry, notably the effective tilt of the fountain 
and the detection inhomogeneity, to reduce the DCP shift 
uncertainty. Based on the experimentally validated model, new 
cavity designs were proposed 1 38 1 to reduce the DCP shift and 
its uncertainty. These measurements lowered the DCP shift 
uncertainty for F02-Cs clock from 3 x 10" ^'^ to 8.4 x 10" 
and defined a set of measurements to establish the DCP 
uncertainty of a fountain, an approach adopted in |40||41|. 
Here, we apply this approach to evaluate the DCP uncertainty 
of FOl, POM and F02-Rb. 

The m = 1 term corresponds to power flowing horizontally 
from one side of the cavity to the other. This contribution can 
be probed by measuring the fountain frequency as a function of 
the tilt of the launch direction. The uncertainty for this term 
is determined by the experimental uncertainty in frequency 
change as a function of tilt, after correctly balancing the feeds 
to null the tilt sensitivity, and by the uncertainty of the fountain 
tilt for nominal operation |39|. When we purposely exaggerate 
the m — 1 DCP shifts by asymmetrically feeding the cavities 
of FOl and FOM, the measurements show tilt sensitivities that 
are much smaller than those of F02-Cs. This is largely due to 
the higher quality factor of their cavities and to their smaller 
cavity apertures. 

Also, when the two feeds are adjusted to produce equal Rabi 
pulse areas, FOl and FOM show no measured tilt sensitivities 
at the level of (1.5 ± 1.7) x lO^^^ mrad^^ and (1.8 ± 1.1) x 
10"^^ mrad~^ respectively, along the axis of the feeds. For 



F02-Rb, this sensitivity is (1.8 ± 1) x 10"^^ mi-ad"i. 

Tilt sensitivities in the perpendicular direction remain to 
be measured for FOl and FOM. Based on the findings of 
||39l , we can reasonably apply, for the time being, the same 
bound as in the parallel direction. For F02-Rb, the measured 
perpendicular sensitivity is < 0.9 x 10"^^ mrad"^. Using the 
same methods as in 1391 (detected atoms versus tilt, frequency 
difference between the 2 asymmetric feedings), we estimate 
the uncertainty on the effective launch direction to be 1 mrad 
for FOl, 0.2 mrad for F02-Rb and 0.7 mrad for FOM in the 
direction of the feeds and 0.7 mrad, 0.7 mrad and 0.7 mrad 
in the direction perpendicular to the feeds. 

The m = and m = 2 terms are calculated from the 
phase distributions Il37lll38l and the corresponding response of 
the fountain. This second step takes as an input a number of 
parameters of the fountain geometry such as the cloud size, the 
cloud temperature and the geometry of the detection. Figure |6] 
shows the measured contrast of the Rabi oscillations at upward 
and downward passages in the microwave cavity, which is used 
to ascertain the fountain geometry to compute the DCP shift 
||39l . The same curve, measured for F02-Cs, is shown in ll42l . 

The resulting DCP shift corrections and their uncertainties 
are presented in Table |l] Several of these uncertainties are 
expected to decrease when improved measurements of the tilt 
sensitivities and tilt uncertainties will be available. 

B. Microwave tensing 

The external degrees of freedom of the atoms evolve dynam- 
ically during the Ramsey interrogation. The interaction of the 
atoms with the standing-wave field in the microwave cavity 
on the upward passage slightly modifies the atomic motion, 
leading to a predicted shift of the clock frequency. This effect 
was first called "the microwave recoil" ||431 . connecting to 
the simple situation where an atom interacts with a plane 
wave as in optical spectroscopy ll44l . In the optical case, 
an atom absorbs a photon from a plane wave. The shift, 
expressed in energy, is (fik)'^ / {2mat), corresponding to the 
kinetic energy of the atom with the momentum fik of the 
absorbed photon. For microwave photons, the corresponding 
fractional frequency shift is 1.5 x 10"^® for ^■^^Cs, comparable 
to the best reported fountain accuracy. 

In 2006 a new approach to this effect provided better 
physical insight and made the calculations tractable to treat 
both transverse dimensions ll45l . The microwave standing 
wave and its spatially varying field amplitude produces a 
radial magnetic dipole force which acts as a lens to focus (or 
defocus) the atomic wave packets, hence the name "microwave 
lensing frequency shift". To calculate the shift with reasonable 
accuracy, it is crucial to properly treat the apertures in the 
fountain. They truncate the downwardly traveling atomic cloud 
and this introduces a difference between the two dressed states, 
which experience opposite lensing deflections on the upward 
cavity traversal. 

Here, we report calculations of the microwave lensing 
for the LNE-SYRTE fountains. We apply the approach that 
was recently used to evaluate the NPL-CsF2 and PTB-CSF2 
fountains BOlBTl. The results are summarized in Table llll The 
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TABLE I 

Corrections and their uncertainties for the distributed cavity phase shift in FOl, F02, and FOM, in units of 10"^^. These 

CORRECTIONS AND UNCERTAINTIES ARE GIVEN SEPARATELY FOR EACH AZIMUTHAL TERM OF THE PHASE DISTRIBUTION, LABELED WITH THE 
AZIMUTHAL number m. The last row is the total DCP correction and uncertainty, GIVEN BY THE COMBINATION OF m = 0, 1 & 2. f FOR 
TECHNICAL REASONS NOT DETAILED HERE, FOl AND F02-Rb ARE TO DATE OPERATING AT A NON-VANISHING PARALLEL EFFECTIVE TILT, HENCE THE 
NON-ZERO m = 1 CORRECTION, f THE FIRST m = 2 CORRECTION FOR FOl IS ZERO WITH A SMALL UNCERTAINTY BECAUSE OF A (45 ± 5)° ANGLE 

BETWEEN THE FEED AXIS AND THE DETECTION AXES. 



simple "analytic" result BOl (P* line in Table [n| neglects all 
apertures on the upward passage and only treats dipole energy 
that is quadratic (cx p^) in the radial coordinate p. Since the 
initial clouds in the LNE-SYRTE fountains are not point-like 
sources, this overestimates the microwave lensing shift because 
it includes atoms that would be cut by the cavity apertures, 
which experience much larger dipole forces away from the 
cavity axis. Including an aperture at the first cavity passage 
eliminates these atoms from the calculation, giving a smaller 
and more realistic shift (2"*^ line in Table ni|i. Adding the 
lowest cutoff aperture for the upward passage, the variation of 
the Rabi tipping angles, and the all-order dipole forces (beyond 
p^), gives a small correction (3'''^ line), as does including 
the detection inhomogeneities and imperfect Ramsey fringe 
contrast, yielding the predicted shift. 

C. Blackbody radiation shift 

The blackbody radiation (BBR) shift is one of the largest 
systematic corrections and a significant source of uncertainty 
in atomic fountains. The bulk of this shift is simply given by 
the scalar Stark shift of the clock transition |461, its response 
to a static electric field. A high accuracy measurement of the 
Stark coefficient was made early on at LNE-SYRTE with FOl 
||47| and subsequently used by most groups for their BBR 
corrections. These measurements where made with applied 
electric fields ranging from 50 to 150 kV.m^^, significantly 
higher than the 831.9 V.m^^ rms value of the BBR field at 
300 K. The greatly improved measurement capability enabled 
by the CSO (see |II-E| i provides an opportunity to measure 
the Stark coefficient at an electric field strength that is much 
closer to a typical BBR field. The early measurements BtII 
were specifically tested for a possible deviation from a purely 
quadratic dependance on the electric field. This test did not 
show any such deviation, which was expected since there 
was no theoretical basis to expect significant higher order 
terms given the measurement accuracy and field strengths. 
Nevertheless, it was important to re-examine these measure- 
ments at much lower fields given the intervening conflicting 



measurements and calculations of ll48lB9l and ll50l . 

These new measurements lISTl at field strengths ranging 
from 1.5 to 25 kV.m^^, are in excellent agreement with 
the earlier measurements of ll47l . To extract the scalar Stark 
coefficient for the BBR clock correction from these mea- 
surements, the small tensor contribution has to be subtracted 
PPTl . In im, a sign error in the earlier theory Il53lll54l 
for this tensor part was found, which was then confirmed 
experimentally |55|. Taking into account this revisited theory, 
we get a corrected value for the scalar Stark coefficient: 
fco = -2.282(4) X 10-1'' Hz.V-2.m2. The change in the BBR 
fractional frequency correction from the previous value of lAlX 
is 6.8 X 10^^*^, 80% of which is due to the sign correction of 
the theory. The uncertainty of fco continues to be dominated 
by the uncertainty on the geometry of the copper plates used 
to apply the electric field in the fountain. The value agrees 
with the most recent theoretical calculations |56||57||38l. It 
also agrees with several other calculations [59|[60|[46|[48J, 
after some Il59lll48l are corrected for omitted contributions 
pointed out in f56l and {5l\. The only discrepancy that seems 
to remain, even after the omitted contributions are included, 
is the calculation of 1521 (see the table summary in [51 1). 
We also note that a thermal atomic beam experiment recently 
provided another Stark shift coefficient measurement |61 1. The 
still preliminary value is in agreement with the LNE-SYRTE 
value and is ~ 5 times less accurate. 

We have also performed direct measurements of the BBR 
shift by inserting a blackbody tube in the FOl fountain [62|, 
exposing atoms to thermal radiation at temperatures ranging 
from 300 to 440 K. Measurements were in excellent agreement 
with the above mentioned scalar Stark shift coefficient, but 
with a 22 times larger uncertainty |5r|. On going comparisons 
between a room temperature and a cryogenic fountain [631 
should yield another measurement of the Cs BBR shift at room 
temperature. 

Regarding the BBR coefficient of the ^^Rb hyperfine fre- 
quency, the experimental status is unchanged from [64 1. On 
the theoretical side, high accuracy ab initio calculations have 
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FOl 


F02-Cs 


FOM 


F02-Rb 


Analytic 


7.8 


8.4 


11.9 


7.0 


Upward cavity aperture 


6.1 


7.0 


9.2 


5.8 


All order dipole & upward cutoff aperture 


5.7 


6.5 


8.6 


5.6 


Detection inhomogeneity & Ramsey contrast 


6.5 


7.3 


9.0 


6.8 



TABLE II 

Calculated microwave lensing fractional frequency shifts for FOl, FOl, and FOM, in units of 10"^''. The microwave lensing 

SHIFT IS GIVEN FOR EACH FOUNTAIN WITH THE APPROXIMATIONS DESCRIBED IN THE TEXT. FROM TOP TO BOTTOM, THE DEGREE OF SOPHISTICATION 

OF THE MODEL IS INCREASING. OUR BEST FINAL RESULT IS GIVEN IN THE LAST ROW. 
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Fig. 6. Measured transition probabilities for resonant Rabi pulses on the 
upward (Rl) and downward (R2) cavity traversals, as a function of the 
microwave amplitude of FOl (top graph) and FOM (bottom graph). 



recently been published in ||65 



D. Microwave leakage 

In the last few years, we have developed a new approach 
to assess frequency shifts due to microwave leakage. Ideally, 
after the state selection process, atoms interact with the 
interrogating microwaves only in the microwave cavity during 
the upwards and downwards traversals. In practice, unintended 
interactions can occur with a residual microwave field that 



is coherent with the probing field. Being unintended, such 
a residual field is by nature uncontrolled, with a potentially 
complex structure with both a large standing-wave and a 
large traveling-wave component. Such an extraneous field, 
which is difficult to exclude at the required level with purely 
electric means, leads to a frequency shift which generally has a 
complicated variation with the fountain parameters, especially 
the microwave power |67 1|68 1. If present, the frequency shift is 
generally difficult to disentangle from other systematic errors 
such as the distributed cavity phase or the microwave lensing 
frequency shifts. 
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Fig. 7. Schematic of a switchable microwave synthesizer with a low phase 
transient using a Mach-Zehnder interferometric switch. 

Our latest experimental approach to evaluate microwave 
leakage shifts uses microwave synthesizers, which not only 
will have intrinsically low microwave leakage but also can be 
switched off when the atoms are outside the microwave cavity. 
The challenge is to develop a switch that does not cause a 
large phase transient of the transmitted field after it is switched 
on. Figure |7] depicts the architecture of our synthesizers. The 
interrogation signal at 9.192 GHz is generated by mixing a 
high power 8.992 GHz with a small amount of 200 MHz signal 
to produce only the minimum needed power. 

This architecture provides an ability to switch the carrier at 
the output by switching the 200 MHz signal. To minimize 
phase transients, the 200 MHz is switched using a Mach- 
Zehnder interferometric switch. The attenuator and the phase 
shifter are adjusted so that, when the PIN diode switch is 
closed, letting the signal pass through, destructive interference 
gives a dark fringe fringe at the output of the interferometer 
When the PIN diode switch is open, the on-state signal passes 
the direct channel of the interferometer with no component 
that is likely to produce a phase transient. Specifically, the 
phase transient occurring in the PIN diode as it is switched is 
circumvented. 

The development of a heterodyne phase transient ana- 
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lyzer ifTOl was crucial to the successful implementation of 
switchable microwave synthesizers in Fig. [7] with low phase 
transients. This analyzer was used extensively to select compo- 
nents and identify the critical parameters to achieve the desired 
performance. Furthermore, this device can also be used to 
unveil spurious phase perturbations that are synchronous with 
the fountain cycle, with sensitivities as low as 1 /irad.s^^ at 
the 9.192 GHz carrier frequency. Such perturbations can occur 
if some elements of the fountain sequence (switching of fairly 
high levels of RP power in acousto-optic modulators, switch- 
ing of the state selection microwave synthesizer, etc) perturb 
the interrogation synthesizer through ground loops or other 
electromagnetic compatibility weaknesses. Required for the 
heterodyne phase transient analyzer, an independent reference 
synthesizer, the 6 — 9 GHz synthesis "for other purposes" in 
Fig. |4] is distributed through an actively compensated fiber 
link. 

With a suitably characterized switchable synthesizer, it 
becomes possible to test for a putative microwave leakage 
by measuring the clock frequency with and without switched 
microwaves [10 1. All the LNE-SYRTE fountains (including 
F02-Rb) are now equipped with a switchable synthesizer. Fur- 
thermore, synchronous perturbations are periodically searched 
for with the heterodyne phase transient analyzer. 



E. Accuracy budget 

Our approach to deal with other systematic shifts have been 
reported elsewhere and have not changed significantly since 
our last reviews Il69lll70l . 

In FOl and F02-Cs, the large cold collision shift |4| is 
extrapolated to zero density in real time through interleaved 
measurements at a high density and a halved density. We use 
the interrupted adiabatic passage method [12J which gives 
accurate ratios of the atomic densities. Owing to the much 
lower atom number dependent effects in ^^Rb 17111172111731 . 
we change the atom number and density with the more 
conventional ways, either loading time or Rabi flopping in the 
state selection cavity. It is also possible to adopt an approach 
where a constant atom number is used most of the time, the 
collision shift being corrected based on a previously mea- 
sured coefficient. The coefficient is measured with dedicated 
comparisons of high and low density, performed a few times 
per year. This approach is also used for FOM, which has a 
comparatively low cold collision shift due to a smaller atom 
number. Note that recently, another method has been proposed 
and demonstrated in order to deal with the cold collision shift 
of Cs I74lll75l . 

The accuracy budget for the LNE-SYRTE fountains is given 



in Table III We give the most recent status, taking into account 



the latest evaluations of the distributed cavity phase shift and 
microwave lensing frequency shifts as reported in Table |l] 
and [n] In the absence of an independent calculation and/or 
experimental investigation, the uncertainty for the microwave 
lensing is conservatively taken to be equal to the calculated 



IV. Applications in time and frequency metrology 
A. ^^'^Cs Fountain comparisons 

For more than a decade, the development of LNE-SYRTE 
fountain ensemble of Fig. |4]has allowed high-accuracy com- 
parisons between atomic fountain clocks. This has especially 
been the case in the last few years during which hundreds of 
days of clock comparisons have accumulated. 

To deal with the differences in clock cycles, mode of oper- 
ation, schedule for verifying the magnetic field, etc, we apply 
the following data processing scheme. Each fountain produces 
measurements, corrected for all systematic biases, at a rate of 1 
measurement per cycle. When an interleaved sequence of high 
and low density measurements is used, this first step requires to 
use these interleaved data to determine the cold collision shift 
coefficient for a particular time window, and then to correct 
the high and low density frequency measurements accordingly. 

The individual fountain measurements are averaged over 
time intervals of 864 s (1/lOOth of a day) starting at the 
beginning of the day. Intervals are discarded if the number 
of valid fountain cycles is below a chosen threshold, for 
instance, due to an automated verification of the magnetic 
field, a laser going out of lock, or a phase/frequency jump 
in the output of the cryogenic oscillator Other filters can be 
applied based on fault of the CSO, one of the fiber links, 
or other subsystems. This leaves a set of validated data where 
the sampling sequence is matched between all fountains. These 
data represent measurements by each fountain of the frequency 
of the common ultra-stable reference, i.e. the 11.98 GHz 
weakly phase locked to one of the hydrogen masers (see Fig. 
|4|, which are used for calibration of the international atomic 
time (see |IV-C| below). 
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value. The last row of Table III gives the status in accuracy 
budgets prior to 2011. 



Fig. 8. Fractional frequency instability between the F02 and FOM Cs 
fountains as a function of the averaging time r. The data corresponds to 
50 effective days of comparisons between MJD 55551 and MJD 55663. Red 
points with en'or bars gives the instability at t = 2^ X 864 s. The black 
curve shows the instability for all r. 

As a third step, the frequency difference between the 
synchronous validated data of two fountains eliminates the 
frequency of the common-mode ultra-stable reference and al- 
lows a direct comparison. An example of frequency instability 
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FOl 


F02-CS 


FOM 


F02-Rb 


Quadratic Zeeman shift 


-1274.5 ± 0.4 


-1915.9 ± 0.3 


-305.6 ± 1.2 


-3465.5 ± 0.7 


Blackbody radiation 


172.6 ± 0.6 


168.0 ± 0.6 


165.6 ± 0.6 


122.8 ± 1.3 


Collisions and cavity pulling 


70.5 ± 1.4 


112.0 ± 1.2 


28.6 ± 5.0 


2.0 ± 2.5 


Distributed Cavity Phase shift 


-1.0 ± 2.7 


-0.9 ± 0.9 


-0.7 ± 1.6 


0.4 ± 1.0 


Spectral purity & leakage 


< 1.0 


< 0.5 


< 4.0 


< 0.5 


Ramsey & Rabi pulling 


< 1.0 


< 0.1 


< 0.1 


< 0.1 


Microwave lensing 


-0.7 ± 0.7 


-0.7 ± 0.7 


-0.9 ±0.9 


-0.7 ±0.7 


Second order Doppler shift 


< 0.1 


< 0.1 


< 0.1 


< 0.1 


Background collisions 


< 0.3 


< 1.0 


< 1.0 


< 1.0 


Total 


-1033.1 ± 3.5 


-1637.5 ± 2.1 


-113.0 ±6.9 


-3341.0 ±3.3 


Prior to 2011 (see note 176,1) 


-1031.4 ± 4.1 


-1635.9 ± 3.8 


-111.4±8.1 


-3340.7 ± 4.2 



TABLE III 

Systematic fractional frequency corrections for FOl, F02 and FOM, in units of 10~^^. The table gives the most recent values 

TAKING INTO ACCOUNT THE RESULTS OF TABLEs[I]AND[n] THE LAST ROW GIVES THE OVERALL CORRECTIONS AND UNCERTAINTIES PRIOR TO 20 1 1 

(details in note f76l ). 



between FOM and F02 over 50 effective days of averaging 
is shown in Fig. l8] Here, the overall short-term instability is 
1.3 X and the measurement resolution is ^ 2 

part in 10^^ after averaging for 4.5 days. The best reported 
instability for this type of comparison is 5.0 x IO^^^t^^/^ 
Il69lll70lll9l. 

Fountain comparisons following the above data processing 
scheme are now performed regularly at LNE-SYRTE. It was 
recently incorporated into a fully automated procedure. In 
addition to daily updates of the fountain comparisons, an 
overview of the systems' parameters is generated to help assess 
the status of the clock ensemble in quasi-realtime with hourly 
updates. To produce fully validated data for fountain compar- 
isons, for example for timescale steering and TAI calibrations, 
the entire data processing is critically scrutinized and redone 
if necessary. Figure [9] shows the monthly average of fractional 
frequency differences D between the three Cs fountain pairs 
from November 2006 to April 201 1. These data coiTespond to 
a total of ^ 1000 effective days of fountain comparison. The 
eiTor bars are the total uncertainties for a given clock pair, 
which are dominated by the systematic uncertainties of the 
two fountains. The overall fractional frequency differences D 
(weighted averages) are reported in Table |IV| together with 
the corresponding overall uncertainties U. \D/U\ represents 
the deviation from a null difference scaled to the la total 
uncertainty. The reduced x^, either with respect to the mean 
(xi)) or with respect to zero (xg) gives a further consistency 
check of these data. Finally, the goodness-of-fij^ as defined 
in ITTTI can be computed, either with respect to the mean 
{Qd) or with respect to zero (Qo)- The Q's very close to 1 
for F0M-F02 and FOM-FOl indicate that fluctuations about 
the means are smaller than the total uncertainties. Qo which 
tests deviations (and/or fluctuations) from the expected null 
difference shows that F0M-F02 and FOM-FOl agree within 
the total uncertainties. Instead the low and Qq for the 



'^The goodness-of-fit Q is the probability that a value of chi-square as poor 
as the value found with the data should occur by chance. Quoting 1771 : "If 
Q is larger than, say, 0. 1, then the goodness-of-fit is believable". 



more stringent F01-F02 comparison would indicate a slightly 
underestimated systematic effect, of a few 10^^®, which is 
under investigation. Note that for the data of Fig. ^\ the 



modifications from our recent study of the DCP shift (see III-A 
and 1 39 1) were not yet implemented. Consequently, the DCP 
shift is one potential candidate to explain these deviations. 
Future comparisons, after a full implementation to all fountains 
of the modified approach to the DCP shift, will clarify this 
point. 

Note that a somewhat similar analysis of worldwide com- 
parisons of fountain primary frequency standards exploiting 
their calibrations of TAI is reported in ifTSl . 
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Fig. 9. Fractional frequency differences between the three pairs of LNE- 
SYRTE Cs fountains from November 2006 to April 20 1 1 . Error bars are the 
total uncertainties which are dominated by systematic uncertainties. 



B. A secondary representation of the SI second 

Following the early demonstration of the low cold collision 
shift in ^''Rb CDCIl, LNE-SYRTE has performed a series of 



11 





Duration (days) 


D 


U 


\D/U\ 


xi 


Qd 


xt 


Qo 


F0M-F02 


320 


-3.0 


9.1 


0.33 


0.28 


0.996 


0.37 


0.983 


FOM-FOl 


240 


4.1 


9.2 


0.44 


0.12 


0.999 


0.28 


0.989 


F01-F02 


500 


-5.7 


6.2 


0.92 


1.48 


0.55 


2.12 


0.001 



TABLE IV 

Summary of Cs fountain comparisons over the November 2006 to April 2011 period. The first column gives the cumulated number 

OF DAYS of comparison. D IS THE FRACTIONAL FREQUENCY DIFFERENCE AND U THE OVERALL UNCERTAINTY, BOTH IN UNITS OF lO"'^®. xi) IS THE 
REDUCED WITH RESPECT TO THE MEAN DIFFERENCE D . Xq IS THE REDUCED X^ WITH RESPECT TO THE EXPECTED NULL DIFFERENCE. Qo AND Qo 

ARE THE CORRESPONDING GOODNESS-OF-FIT AS DEFINED IN 



high accuracy absolute measurements of the *'^Rb hyperfine 
frequency against Cs fountain clocks |[64lll79l ll3]|. The 2003 
comparison was the most accurate of any atomic frequency 
comparison at the time. Since then, it has been superseded 
by the rapidly improving optical frequency measurements. 
The ^^Rb hyperfine frequency was the first secondary rep- 
resentation of the SI second recognized by BIPM |80|. Its 
recommended value is given by the 2002 LNE-SYRTE mea- 
surement: 6 834 682 610.904324 Hz with a recommended 
uncertainty of 21 /xHz or 3 parts in 10^^. This recommended 
uncertainty was chosen to be ^ 3 times larger than the 
actual uncertainty of the measurement (1.3 x 10~^^). Our 
latest reported measurement (4 /iHz uncertainty or 5.9 parts 
in 10^^), together with an analysis of the consistency of all 
our measurements can be found in |3|. An independent ^^Rb 
versus ^'^^Cs fountain comparison has been reported for the 
first time in 2010. Agreement with LNE-SYRTE measure- 
ments is claimed at the 10~^^ level but the actual result of 
the measurement has not yet been reported iSTI ll82l 183 1 . Also 
notable is the construction of at least 8 Rb fountains at the 
US Naval Observatory ||84|li85l . 



C. Contributions to TAI 

Since 2007, LNE-SYRTE fountains have measured the fre- 
quency of one of LNE-SYRTE H-masers almost continuously. 
These data are used to calibrate the International Atomic 



mercial Cs clocks located in France and elaborated at LNE- 
SYRTE. This steering has improved the long-term stability of 
the timescale with respect to TAI ifSTl . 



2005 2006 2007 2008 2009 2010 2011 



Time (TAI). Figure 10 shows all calibrations of TAI by LNE- 
SYRTE cesium fountains from November 2004. Each data 
point corresponds to a measurement, by one fountain, of the H- 
maser average frequency over a 20—30 days period, reaching a 
typical statistical uncertainty of 3x 10^^^, which is augmented 
in quadrature by an uncertainty due to dead times of typically 
1.5 X 10-16 ill. Between the 2007 and 2010 period, LNE- 
SYRTE reported 78 formal evaluations to BIPM, out of a 
total of 161 fountain reports worldwide, i.e. LNE-SYRTE 
has contributed almost 50% of all reports since 2007. For 
contributions to TAI, the gravitational redshift corrections are 
-69.3 X 10-l^ -65.4 x lO'^^ and -68.7 x lO-^^ for FOl, 
F02 and FOM respectively, with an uncertainty of 1.0 x 10^ 
corresponding to a 1 m uncertainty of height above the geoid. 
Contributions to TAI also provide a vehicle to compare to 
other fountain clocks around the world. The latest update on 
these comparisons appears in 1781 . 

Since mid-2006, fountain data, including from F02-Rb, are 
used to steer the French Atomic Time TA(F) on a monthly 
basis. TA(F) is a timescale based on an ensemble of com- 




53200 53600 54000 54400 54800 55200 55600 
MJD 

Fig. 10. Calibrations of the International Atomic Time (TAI) by the LNE- 
SYRTE fountains since 2004. The data are froin BIPM circular T. Each point 
is the fractional frequency difference between TAI and the SI second as 
measured by one of LNE-SYRTE's fountains and corresponds to a formal 
monthly report to BIPM. The red curve is the difference between TAI and 
the SI second, as computed by BIPM from all available primary standards. 



D. Absolute measurements of optical frequencies 

The accuracy and stability of optical atomic clocks have 
now largely surpassed those of microwave clocks, and notably 
atomic fountain clocks |88||89||90||91 1. This opens a path to 
redefining of the SI second with an optical transition |92|. 
One of the prerequisites to a redefinition is absolute frequency 
measurements of potential optical transitions with the highest 
accuracy. The LNE-SYRTE fountain clock ensemble (Fig. |4]) 
offers some unique features to contribute to this task. This 
includes the availability of several fountains for redundant 
measurements, of an ultra-low noise reference signal that pro- 
vides the best reported stability to date between a microwave 
and optical clock |93|, and of the transportable fountain FOM 
for measuring optical frequency at other laboratories. Abso- 
lute frequency measurements of optical frequencies can also 
probe the stabiHty of fundamental constants Ii94lll95lll96lll97l . 
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determine the Rydberg constant 11981 and test Quantum Elec- 
trodynamics (QED) |99l. 

Next, we give an overview of absolute optical frequency 
measurements made with LNE-SYRTE fountain ensemble. 

1) Hydrogen IS — 2S transition measured with FOM at 
MPQ Garching, Germany: At the Max Planck Institut fiir 
Quantenoptik (MPQ) in Garching (Germany), we have mea- 
sured the absolute frequency of the 1S-2S optical transition 
via 2 photon excitation at 246 nm |100||95||1()1J. Three mea- 
surement campaigns were completed in 1999, 2003 and 2010. 
For these experiments a femto-comb laser is used to measure 
the frequency of the excitation laser. The repetition rate is 
locked to a hydrogen maser whose frequency is simultaneously 
measured directly by FOM. The noise of the maser signal is 
filtered by phase-locking a BVA quartz oscillator to improve 
the measurement stability. The last measurement in 2010 had 
the best fractional uncertainty: 4.2 x 10~^^ 1101 1. The mea- 
sured 1S-2S transition frequency is 2 466 061 413 187 035(10) 
Hz. 

High resolution spectroscopy of hydrogen is interesting to 
test highly accurate atomic structure calculations for QED 
tests M102I . Also, the recent capture of antihydrogen at CERN 
01O3I1 is a first step toward cooling of antimatter to enable 
high resolution spectroscopy of antihydrogen for fundamental 
tests of charge, parity, time reversal (CPT) symmetry and the 
gravity of antimatter. 

2) ^"Ca^ optical clock measured with FOM at the Uni- 
versity of Innsbruck, Austria: In 2007 at the Institut fiir 
Experimentalphysik and Quantenoptik, University of Inns- 
bruck (Austria), we measured the absolute frequency of the 
4s ^5*1/2 — 3d ^-05/2 optical quadrupole transition of a 
trapped single *°Ca+ ion. The frequency of this transition is 
411 042 129 776 393.2 (1.0) Hz fEl- This coiTesponds to a 
fractional frequency uncertainty of 2.4 parts in 10^^. During 
this experiment, FOM was used in the autonomous config- 



uration mentioned in II-D where a BVA quartz oscillator is 
frequency locked to the FOM spectroscopic signal. The output 
of the BVA quartz oscillator synchronizes the repetition rate of 
the femto-comb laser that measures the probe laser frequency 
for the ^°Ca+ transition. This quadrupole transition in ^°Ca+ 
ion is another candidate for an optical clock and for tests of 
the stability of fundamental constants. 

3) ^^Sr and *^Sr optical lattice clocks measured at 
SYRTE.- LNE-SYRTE is developing two Sr optical lat- 
tice clocks. During their development, several absolute fre- 
quency measurements of the Sr clock frequency were made 
ll93]llfb5]|||T06llllT07ll . The best reported uncertainty for these 
measurements is 2.6 x 10~^^ ||93|. Both Titanium: Sapphire 
based and Er-doped fiber based optical frequency combs 
were used to measure the frequency of the 698 nm laser 
light stabilized to the Sr atoms. The 2006 measurements of 
the ^^Sr used a fiber optical frequency comb transported 
from the Physikalisch Technische Bundesanstalt (PTB) in 
Braunschweig (Germany) |106||93|. As already mentioned, 
the ultra-stable reference based on a CSO (Fig. |4]l allows 
the best reported short-term stability to date for comparing 
an optical clock to a microwave clock: 6 x 10~^^r^^/^. 
In the recent measurements, the optical frequency comb is 



stabilized to an ultra-stable laser as described in ll22l . making 
the repetition rate of the femtosecond laser an ultra-stable 
microwave signal. This signal is then measured against the 
CSO based ultra-stable reference which is in turn measured 
by the fountain clocks. Simultaneously, the beat between the 
optical frequency comb and the 698 nm light is counted to 
extract its absolute frequency. Combining these with similar 
measurements by other institutes, stringently tests the stability 
of fundamental constants |94|. 

4) Measurement of ^^^Hg and ^°^Hg optical clock tran- 
sition frequencies at SYRTE.- LNE-SYRTE is developing 
an optical lattice clock based on neutral mercury. The first 
absolute frequency measurements of the ^^^Hg and ^°^Hg 
optical clock transition were in 2008 |1081 . These initial mea- 
surements improved the previous knowledge of the transition 
frequency by more than 4 orders of magnitude, to a fractional 
frequency uncertainty of 5 x 10~^^. More recently, measure- 
ments against LNE-SYRTE ultra-stable reference were used 
to perform the first experimental determination of the magic 



wavelength for Hg II109I . Ongoing development of the Hg 
clock is expected to soon yield measurements at the 10~^^ 
level and beyond. 

E. Other measurements with the transportable fountain FOM 
We have performed frequency comparisons with the engi- 



neering model of the space clock PHARAO (see VI below and 
1 1101) between 2007 to 2009, at the Centre National d'Etudes 
Spatiales (CNES) in Toulouse, France. We have also verified 
the frequency performance of the ACES architecture pi 1 \\. 
The test included a ground model of the onboard data handling 
unit (XPLC), the Frequency Comparison and Distribution 
Package (FCDP), the PHARAO engineering model, and a 
ground model of the Space Hydrogen Maser (SHM). The two 
clocks are combined to generate a timescale with the short- 
term stability of SHM and the long-term stability and accuracy 
of PHARAO. 

F. Remote comparisons via T2L2 satellite laser link 

In 2010, time transfer by the T2L2 laser link via the 
JAS0N2 satellite 111211 was tested by comparing FOM at the 
Observatoire de la Cote d'Azur (OCA) in Grasse (France) 
with the other LNE-SYRTE fountains at the Observatoire de 
Paris. The common time transfer techniques, carrier phase 
GPS and TWSTFT, were also used in parallel. Notable in this 
comparison is the large gravitational redshift, — 1.384 x 10"^'^, 
at the OCA altitude of 1268 m. Preliminary results |113| show 
that all the time transfer methods are consistent to within 2 ns 
over 2 months. The frequency difference between the remote 
fountains is measured with a typical uncertainty of 1 x 10~^^, 
which determines the gravitational redshift difference to 10~^. 
It now remains to refine the analysis and to evaluate the 
ultimate performances of the T2L2 time transfer link. Note that 
at MPQ, at CNES and at OCA, FOM also provided calibration 
of the TAI by using a GPS receiver and a carrier phase GPS 
analysis software developed by CNES or by NRCAN. This is 
the first time that a primary frequency standard contributes to 
steering the TAI from different sites with different gravitational 
redshift II114I. 



13 



V. Applications in fundamental physics 

One of the most interesting applications of atomic clocks 
is to test fundamental physics. In this section, we review 
contributions of LNE-SYRTE fountain ensemble to such tests. 



A. Stability of constants 

Repeated high accuracy comparisons between atomic (or 
molecular) frequencies can test the stability of fundamental 
constants such as the fine structure constant a or the electron 
to proton mass ratio rrie/mp (see, for instance, f 1 151 Pl 161 ™d 
references therein). Laboratory experiments usefully comple- 
ment tests over cosmological timescales since the interpreta- 
tion does not rely on any assumptions about a cosmological 
model JTlTllfTTSl . 

A first test performed with LNE-SYRTE fountain ensemble 
comes from a series of ^^Rb versus ^■^•^Cs hyperfine frequency 
comparisons made during the development of the F02 dual 
Rb/Cs fountain |79||69||7()1 and llTSl. In this last 2008 report, 
a putative time variation of the ratio of the two hyperfine 
frequencies is constrained to (iln(z/Rb/i'Cs)/'^^ — (—3.2 ± 
2.3) X 10~^^ yr~^. In terms of fundamental constants, this 
result yields d\n{a-°-'^^[gYib/gcs])/dt = (-3.2 ± 2.3) x 



LNE-SYRTE fountain's contributions to this test were the 



10~ yr~ where g^b and gcs are the nuclear g-factors 
in ^^Rb and ^'^'^Cs. Expressing these g-factors in terms of 
fundamental parameters of the Standard Model M120H121II 
gives: dln(a-0-49[m^/AQCD]-° "2^)/di = (-3.2 ± 2.3) x 
10~i^ yr~i, where Aqcd is the mass scale of Quantum 
Chromodynamics (QCD). An improved analysis including 
more recent measurements, with F02 in the dual fountain 
configuration, will improve this value significantly. 

A second test comes from a series of absolute fre- 
quency measurements of the ^^Sr optical lattice clock 
by LNE-SYRTE (see |IV-D3| l, the University of Tokyo, 
Japan and JILA, Boulder, Colorado, USA [94 1. To- 
gether, these constrain the putative variation of I'si/i^Cs 
to dh-i{vsr/vcs)/dt = (-7 ± 18) x IQ-^^ yr^^ coitc- 



sponding to dln{a 



2.77 r 



j/Aqcd] 



-0.039 



me/AQCB])/dt = 



(7 ± 18) X 10^^° yr^ . The same measurements con- 
strain the putative variation of this same combina- 
tion of constants with the gravitational potential to 
c2dln(a2.77[^^/AQCD]-°°^9[m,/AQCD])/rfC/ = (-5.8 ± 
8.9) X 10~^, where c is the speed of light and U the grav- 
itational potential |^ This test relies on exploiting the yearly 
modulation of the gravitational potential of the Sun due the 
eccentricity of the Earth orbit. 

Measurements of the H(1S-2S) transition performed at MPQ 



Garching (see IV-Dli using the transportable fountain FOM 



as a reference, offer a third test: (iln(fc'H(is-2S)/'^Cs)/rf^ = 
(-3.2 ± 6.3) X 10-^"" yr-i |95|. This translates to 

dln(Q2.«3[^^/AQCD]-° °3^K/AQCD])M - (3.2 ±6.3) x 
10~i5 yr^i, a constraint which will soon improve as a result 
of a recent measurement campaign |101|. 

A fourth test comes from a series of comparisons be- 
tween the ^^•^Cs and the hydrogen hyperfine frequencies. 

^For instance, the gravitational potential created by a point mass m at 
a distance r is: U{r) = Gm/r where G is the Newtonian constant of 
gravitation. 



calibrations of TAI reported in BIPM Circular T (see IV-C I 
enabling the connection with 4 hydrogen masers participating 
to the elaboration of the NIST ATI atomic timescale (see 
for instance fT22|). Other Cs fountain contributions to the 
experiment came from PTE in Germany and from the Istituto 
Nazionale di Ricerca Metrologica (INRIM) in Italy. The anal- 
ysis was performed at the National Institute of Standards and 
Technology (NIST) Boulder, Colorado, USA 1 123 1 and gives 
\c^d\n{vii/h'cs)/dU\ = (0.1 ± 1.4) x 10"'^, corresponding to 
|c2dln(aO-83[r7i,/AQCD]°-")/rfC/| = (0.1 ± 1.4) x IQ-^. 

Note that instead of using a, mg/AqcD and TOe/Aqcd 
as independent variables, a, fj, — me/nip and mq/irip 
can be used. The link between the two approaches is 
dln(77ip/AQCD) ^ 0.048 x dln(m,/AQCD) LHU- 



B. Test of Lorentz Invariance in the matter sector 

The LNE-SYRTE F02 fountain has tested the anisotropy 
of space as in Hughes-Drever experiments (see for instance 
111241 ). For this experiment, the fountain sequence is tailored 
to probe opposing Zeeman transitions in ^■^•^Cs to test for a 
putative variation of their frequencies when the orientation 
of the quantization axis changes (here, due the Earth rotation) 
with respect to a supposedly preferred frame [125], such as the 
frame of the Cosmic Microwave Background. The experiment 
was interpreted within the framework of the Lorentz violating 
Standard Model Extension (SME), where it is sensitive to pro- 
ton parameters corresponding to a largely unexplored region of 
the SME parameter space. The constraints for 4 parameters, 
already constrained by other measurements, were improved 
by as much as 13 orders of magnitude and 4 parameters were 
constrained for the first time. Operating F02 in the dual Rb/Cs 
configuration f3l opens new possibilities for improving and 
complementing these tests. 

C. Other tests 

The development of LNE-SYRTE fountain ensemble also 
offered a test of Lorentz Local Invariance in the photon 
sector, not using the fountains, but simply comparing the CSO 
to a hydrogen maser (see Fig. |4]l for a duration that now 
exceeds 10 years. This experiment is to date the most stringent 
Kennedy-Thorndike test by a factor of 500. The latest update 
is published in II126I . The experiment can also be interpreted 
in the SME framework lfT27]l . 

VI. Development of the PHARAO cold atom space 

CLOCK 

Atomic Clocks Ensemble in Space (ACES) is a European 
Space Agency (ESA) fundamental physics mission originally 
proposed by Laboratoire Kastler Brossel and SYRTE. It is 
based on the operation of highly stable and accurate atomic 
clocks in the microgravity environment of the International 



Space Station (ISS) 111 1 II . The time scale generated by the 
ACES clocks on-board the ISS is delivered to Earth through 
a high-performance two-way time and frequency transfer 
link. The clock signal is used to perform space-to-ground as 
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tested the engineering model (PHARAO-EM) of the clock 



Fig. 11. The engineering model of the PHARAO cold atom space clock 
under test at CNES, Toulouse (photo is courtesy of CNES). The picture shows 
the vertically mounted "cesium tube", the vacuum chamber surrounded by 
the magnetic shields. The gold-coated box behind the cesium tube contains 
the "laser source". A second gold-coated box, below the laser source, is the 
microwave synthesizer. The black box in front, at the bottom of the cesium 
tube, is the on-board computer unit. 



well as ground-to-ground comparisons of atomic frequency 
standards. The ACES scientific objectives cover both funda- 
mental physics and applications. Tests of Special and Gen- 
eral Relativity will be performed with improved accuracy, 
as well as a search for temporal variations of fundamental 
constants. On the application side, frequency comparisons 
between distant clocks, both space-to-ground and ground- 
to-ground, will be performed worldwide with unprecedented 
resolution. ACES will demonstrate a new type of "relativistic 
geodesy", which, based on a precision measurement of the 
Einstein's gravitational red-shift, will resolve differences in 
the Earth gravitational potential at the 10 cm level. Finally, 
ACES will contribute to the improvement of global navigation 
satellite systems (GNSS) and to future evolutions of these 
systems. It will demonstrate new methods to monitor the 
ocean surface based on scatterometric measurements of the 
GNSS signal and it will contribute to the monitoring of the 
Earth atmosphere through radio-occultation experiments. The 
expected performance is a time stability of 10 ps over ten 
days (1 ps for common view comparisons between ground 
clocks) and a frequency accuracy better than 3 x 10~^^. The 
accuracy and the long term frequency stability are defined by 
the PHARAO clock, an instrument developed by the French 
space agency CNES and SYRTE. The clock uses cold cesium 
atoms, slowly moving through a Ramsey cavity. We have fully 



111 101 . Figure 1 1 shows a picture of the PHARAO-EM under 
test at the CNES assembly room in Toulouse. For nominal 
operation the frequency stability is 3.3 x 10~^^t~^/'^. The 
central Ramsey resonance linewidth is 5.6 Hz and 10^ atoms 



are detected. The contrast is only 93% (Fig. 12 1 since, on 
Earth, the atoms do not spend the same time in the two 
interaction zones of the Ramsey cavity. 
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Fig. 12. Experimental (red) and calculated (blue) Ramsey fringes of the 
PHARAO-EM tested on Earth. The central fringe does not have the highest 
contrast because, operating on Earth, gravity leads to different interaction 
times on the two cavity interactions. 

On Earth the clock cycle duration lasts 790 ms with a 
microwave interrogation time of 90 ms. Consequently the 
phase noise of the microwave signal, synthesized from a 
quartz oscillator, contributes to the frequency instability at a 
level of 2.2 x lO^^'^r^^/^. We have performed an accuracy 
evaluation of the clock. The main frequency corrections and 
their uncertainties are shown in table [V] The total frequency 
uncertainty is 1.9 x 10^^^. The largest uncertainty is the 
quadratic Zeeman shift because the magnetic field along the 
atomic trajectories is not homogeneous due to a problem with 
the innermost magnetic shield. The second largest contribution 
is the cold collision frequency shift. The frequency shift per 
detected atom is similar to other cold cesium clocks and the 
uncertainty is mainly limited by the measurement duration 
(only two weeks). The blackbody radiation shift is deduced 
from measurements of the ambient temperature. A technical 
problem prevented us from using the calibrated temperature 
probes located next to the Ramsey cavity. Finally, the first 
order Doppler frequency shift has a large value since the 
atom's axial velocity decreases due to gravity field, which 
causes the longitudinal DCP shifts from the 2 interactions 
to not cancel. This gives a large DCP shift that has a large 
dependance on the microwave field amplitude. 

The PHARAO-EM clock, operated in several configura- 
tions, was compared to the transportable fountain FOM. After 
applying all systematic corrections (Table III and |V]i, the 
frequency differences were smaller than 2 x 10^^^, with a 





Correction 


Uncertainty 


Quadratic Zeeman effect 


-5633 


14 


Blackbody radiation 


162.5 


2.2 


Cold collisions & cavity pulling 


63 


11 


First order Doppler 


-33 


5 


Total 


-5440 


19 



TABLE V 

Systematic fractional frequency corrections for the 
pharao engineering model during the 2009 tests, in units of 

io-i«. 



resolution of 5 x 10^ . During these measurements, the 
accuracies of the clocks were 8.1 x 10^^^ for FOM and 
1.9 X 10"^^ for PHARAO-EM. Consequently, there is no 
significant frequency offset between the clocks. Testing the 
engineering model demonstrated the full operation of the 
PHARAO clock architecture and highlighted some improve- 
ments to implement in the flight model. The PHARAO flight 
model is now being assembled and it will be delivered in 
2012 together with the other ACES instruments. The launch 
date will be 2 years after the delivery of all the instruments. 

SYRTE will be a master station to survey and to analyze 
the ACES signal. One of the microwave link station developed 
for ACES will be installed on the roof of a building at 
the Observatoire de Paris. Hydrogen masers steered by the 
atomic fountains and UTC(OP) will synchronize the station, 
forming the ACES ground segment at SYRTE. Strontium and 
mercury optical lattice clocks at SYRTE 1 128ll T09l . which 
are connected to the fountain ensemble via optical frequency 
combs, will also contribute to the ACES mission. Notably, 
ACES will provide ground-to-ground comparisons of these 
optical clocks to similar clocks at other institutes around the 
world, demonstrating relativistic geodesy. 

VII. Conclusions 

This paper gives an extensive overview of the develop- 
ment of the LNE-SYRTE fountain ensemble. We can also 
identify some priorities for the future. For fountains, it will 
be important to further improve inter-comparisons to reduce 
inaccuracies to 1 to 2 parts in 10^®, and to possibly re- 
move yet unidentified sources of uncertainty or instability. 
Implementing more sophisticated microwave cavity designs is 
one interesting approach. Simultaneously, the reliability of the 
fountain ensemble must further improve. This is required by 
the ground segment of the ACES mission - quasi-continuous 
operation at the best accuracy of the three fountains for at least 
18 months. This will be achieved by a number of technical 
improvements in the fountains (notably the laser systems), the 
fountain environment, and by complementing the cryogenic 
sapphire oscillator with a second low-noise microwave source 
derived from an ultra-stable laser This system will also ensure 
a permanent link to the optical clocks. Achieving these goals 
will enable new or improved applications: realization of an 
ultra-stable local timescale, continuously steered to the atomic 
fountains, even more regular calibrations of TAI, improved 
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measurements of optical frequencies, enhanced tests of fun- 
damental physics, measurement of the gravitational red shift 
within the ACES project, and remote clock comparisons via 
ACES. 
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